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1 Introduction 



At LEP 2, two-photon collisions make up by far the largest class of events. These are processes 
in which the incoming electrons each radiate a photon, which collide to produce a hadronic 
or leptonic final state. A photon can obviously interact electromagnetically with any charged 
object. However, since it has the same quantum numbers as a vector meson, it can fluctuate into 
one, and can therefore also be considered as an incoming hadron, interacting strongly through 
its partonic constituents. The interplay between these two ways of interacting is unique to the 
photon and provides much of the interest in 77 physics. 

The two-photon invariant mass spectrum is peaked at low mass, so the bulk of events 
only produce a few particles, and resonances and exclusive final states can be studied. The 
total cross-section is so large however, that there are enough events to study deep inelastic e7 
scattering and high-p^ jets and heavy quark production in 77 collisions. 

LEP 2 is the highest energy and luminosity e7 and 77 collider available, and many of the 
same studies as at electron-hadron and hadron-hadron colliders can be made here. However, 
one essential difference is that, since the beam remnants (the electrons and positrons) typically 
leave the detectors undetected, the energies of the incoming photons are not known and must 
be reconstructed from the properties of the final state. In events in which most of the final-state 
particles are visible in the detector, this is easily done. However, at high energies the final-state 
distribution becomes increasingly forward-peaked and much of the energy goes into the very 
forward parts of the detector, or is even missed in the beam pipe. It is therefore essential that 
as much of the final state as possible is measured, in particular to detect hadrons in the forward 
detectors that have hitherto only been used for electron tagging. 

It is also essential that we are able to understand the details of the multi-particle final 
states in these interactions, which puts very high demands on the event generators used in 
the analysis. During the course of this workshop, some of the standard general purpose event 
generator programs on the market have been developed to also handle e7 and 77 collisions. 
This means that we can use our experience from ep and pp collisions to give a more complete 
description of the hadronic final state of two-photon collisions. These models can now be tested 
at LEP 1 and should give reliable extrapolations to LEP 2 energy. 

In section ^ of this report, we will describe briefly the models of the generators we have 
studied during the workshop. Then, in section § some comparisons between the programs are 
presented for different classes of events. In section ^, the main programs are presented in some 
detail, followed by section |^ where some other generators are presented more briefly. Finally 
in section ^ we present our conclusions. Related work is presented in the reports of the "77 
Physics" [H and "QCD Event Generator" p| working groups. 
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2 General Features 



The event generators used for 77 physics can be divided into two main groups. One deals with 
low multiplicity final states, like resonances, exclusive channels and leptonic channels and the 
other with high mass multi-particle final states. 

When there are few particles in the final state, the fully differential cross-section for a given 
process can usually be derived directly from a model of that process. Event generators can then 
be viewed as a particularly convenient numerical implementation of that cross- sect ion, in which 
arbitrarily complicated phase-space cuts and detector simulation can be incorporated. In this 
group are 'four-fermion generators', which incorporate the full set of QED matrix elements for 
e"'"e~ — >■ e"'"e~ff, and more general programs that use a 77 luminosity function to separate the 
process e"'"e~ — > e^e~X into two stages, e^e~ — > e^e~77 and 77 —>■ X. 

The other group of generators describe multi-particle production, for which cross-sections 
cannot be directly calculated in quantum mechanics. They use semi-classical probabilistic 
models to separate the process into several phases. First photons are radiated from the incoming 
electrons to give beams of quasi-real photons. Then a hard sub-process is generated, using 
partonic 2 — > 2 matrix elements folded with the parton densities of the photon. The emission 
of additional partons from the incoming partons is generated by evolving them "backwards" 
in an initial-state parton shower, and from outgoing partons by generating a final-state parton 
shower. Finally, the partons are converted to hadrons, which are then allowed to decay. 

2.1 Photon generation 

A two-photon reaction can be factorized into photon fiuxes of radiation from incoming e^ and 
the final state of two-photon collisions. The decomposed differential cross-section for 

e+(pi)e~(p2) e+(p'i)e"(p2)7(gi)7(g2) ^ e+(p;)e"(p2)X(gi + gs) (1) 

da = f /^^^f,' ~ f\ Y' Upt^pt^crrr + 2|p+-pMr^^ cos(20) 

IGTT^qiqi \{piP2) — 'mfmi / 

+^Pl P2'^rL + 2piP2 (^LT + Pi P2 (^LL - o\Pi P2 rTLCOS[(p)j— —, [Z) 

where the a's and r's are linear combinations of the cross-sections for 77 X of transverse(T) 
and longitudinal (L) photons, the flux factor p"'' has photon helicities labelled by -|-, — , 0. Some 
dedicated 77 generators such as TwoGAM use the full form of Eq. |^, while most models 
simplify further by taking the approximation. At ^ 0, the photons are quasi- real 

and transversely polarized and after integration over 0, the angle between lepton scattering 
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planes, the only remaining term is arr- Expressed in terms of a luminosity function we have 



where Wi = Eb — E[ is the photon energy in the lab frame. The luminosity function can be 
decomposed as the product of a factor for each of the photons, 

(PC (fC 



dwidu2d9id02 duoidOi duj2d92 

where, apart from a trivial kinematic factor, d^C is the Equivalent Photon Approximation 
(EPA) flux factor, 

with X the light-cone momentum fraction and the photon virtuality, = |g^|. As discussed 
in Ref . , there are two important corrections to the usual EPA formula. The first is to include 
the sub- leading term of relative order ml/P^, 

which can give corrections of order 10% for untagged and anti-tagged cross-sections. At present, 
of the QCD event generators, only Phojet includes this correction. The second important 
correction is to include the correct, pro cess- dependent, dynamic upper limit on . However, 
this is only important for untagged cross-sections, and when an anti-tag condition is imposed, 
as it is throughout this report, this corrections become small. 



2.2 Photon distribution functions 

In resolved-photon processes we need parametrizations of the distribution functions for partons 
inside the photon. These obey an inhomogeneous form of the usual evolution equations. As 
discussed in more detail in the report of the "77 Physics" working group f^, their solution 
can be written as the sum of a hadronic or VMD part, which evolves according the usual 
homogeneous equation, and a point-like or anomalous part. 

There are a number of parametrizations available for the parton distribution functions of on- 
shell photons. Most of them are contained in a single package, PDFLIB |^ , to which most of the 
event generators are interfaced. At present none use the recent models of the structure of virtual 
photons, although HERWIG does implement a simple P^-suppression model. In the following 
when comparing different generators, we use the fairly similar SaS ID or GRV LO sets. 
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single resolved 


7q - 


qg 




7g - 


qq 


double resolved 


qq' - 

gg - 


qq' 
qq 


DIS 


eq - 


eq 



Table 1: The standard hard sub-processes used by event generators for different event classes 
2.3 Hard sub-processes 

Having defined the structure functions and parton densities, we can now use the same ma- 
chinery as for pp or ep collisions to generate hard sub-processes, with the exception that we 
here have additional processes where the photon couples directly in the hard interaction. In 
77 collisions we therefore talk about three kinds of events, direct, single-resolved and double 
resolved, depending on whether the photons couple directly or not. In Table |1| the standard 
hard sub-processes are listed for different event classes. Some programs, like Pythia, make 
a further distinction between the anomalous and VMD-like part of the resolved photon and 
therefore have six different event classes. 

In deep inelastic e7 scattering, the exchanged photon is usually more virtual than the struck 
quark, so the EPA is no longer a good approximation (in other words the pro cess- dependent 
dynamic upper limit on mentioned in Section is exceeded). One therefore needs, in 
principle, to use the full 2 — >• 3 processes eq eqg, eg — > eqq and e7 eqq. However, when 
the quark line is much less virtual than the photon, it can be approximated by the DGLAP 
probability distribution to find the quark inside a higher-x quark, q qg, gluon, g — qq or 
photon, 7 — i> qq and hence can be absorbed into the evolution of the photon distribution 
functions. Thus we are again left with a 2 — 2 process, eq — >■ eq, for which one uses the 
lowest-order matrix element. 

Other processes that are usually treated separately are the ones involving heavy quarks, 
where the matrix elements are different from the massless case. Here there are also event 
generators that use the next-to-leading order matrix elements, incorporating one additional 
parton. These can be considered as an exact treatment of the first step of a parton shower 
and can be compared with the usual algorithms discussed below, which contain approximate 
treatments of all steps. 



5 



2.4 Soft processes 



The cross-section for quasi-real 77 scattering is dominated by processes in which there is no 
hard scale. Several models exist for the subdivision of the soft cross-section into separate 
components, typically: elastic, 77 VV; single diffractive dissociation, 77 VX; double 
diffractive dissociation, 77 — >• XX; and inelastic, 77 — >• X. The cross-section for each is given 
by the model, and the total cross-section is simply their sum. Clearly the separation is model 
specific, for example the difference between diffractive dissociation and inelastic scattering is the 
presence of a rapidity gap between the two systems in the former and not in the latter and the 
cross-sections depend on the size of this gap, and only the sum of the processes can be directly 
compared between models. Nevertheless, the available models use fairly similar definitions 
and comparisons between components can prove useful. Since most of these reaction types 
cannot be calculated from first principles, they are characterized by a rather large number 
of adjustable parameters. Nevertheless, since the models assume photon-hadron duality and 
hadron universality, parameters can be fixed in hadron-hadron and photon-hadron collisions, 
giving parameter-free predictions for LEP 2. Phojet and Pythia contain rather complete soft 
interaction models, while HERWIG contains only the non-diffractive part of the cross-section. 

2.5 Multiple interactions 

At increasing centre-of-mass energies most sub-process cross-sections grow faster than the total 
cross-section, eventually overtaking it. Important examples include (supercritical) soft pomeron 
exchange and hard two-parton scattering above any given ptmin cut. This corresponds to the 
possibility of multiple soft or hard scatterings within a single 77 collision. At present only Pho- 
JET implements both soft and hard multiple scattering in the same package with ptmin forming 
the boundary between the two. Pythia and HERWIG (through an interface to the Jimmy 
generator) both implement multiple hard scattering above ptmin, but with soft models that do 
not vary with ptmin, so that it becomes a critical parameter of the model. Although the general 
ideas of the models are similar, there are many specific differences in their implementation. For 
example, the current models only allow multiple interactions in the hadronic cross-section, but 
since the definition of anomalous and hadronic events differs (see below), so do the multiple 
scattering results. In all three cases however, the models are essentially identical for 77 and 
7P collisions, so experience gained at HERA will certainly help constrain the predictions for 
LEP 2. 

2.6 Part on showers 

The hard scattering disturbs the colour fields of the incoming partons and as a result they 
partially 'shake off' the cloud of gluons that normally surrounds colour charges. This gives rise 
to a shower of accompanying radiation, which is conventionally modelled as a series of emissions 
from the incoming parton, starting from the parton entering the hard interaction and tracing 
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its history back towards the incoming photon. During this 'backwards evolution', the emission 
at each step is required to be at a lower scale than the previous one until, near the incoming 
photon, no more radiation is resolvable above the infrared cutoff. This procedure is governed 
by the splitting functions of the DGLAP evolution equation and is guided by the input set of 
parton distribution functions. 

This is the way it is done in the two standard generators HERWIG and Pythia where the 
main difference between the two is the choice of evolution, or ordering variable - Pythia uses 
the virtuality of the evolving parton, while HERWIG uses a generalized emission angle, which 
incorporates colour-coherence effects. 

Another important difference between the two programs is the way they distinguish between 
the anomalous and VMD-like part of the structure functions. In Pythia this distinction is 
done beforehand, and the shower is identical in both cases, only using different parton densities 
corresponding to the two parts of the structure functions. In HERWIG, however, the two parts 
are treated together, and an additional branching is introduced into the parton shower, where 
a quark may be evolved back to the incoming photon. If this happens before the shower is cut 
off, the event is called anomalous, if not, it is a VMD-like event. In both cases, the photon 
remnant will have a larger transverse momentum in the anomalous case — this is automatic in 
the perturbative treatment of HERWIG, while in Pythia it is put in by hand. 

In both programs the generated partons may continue to branch in a final-state shower, in 
the same way as is done in e"'"e~ annihilation. This is discussed in more detail in the report 
from the "QCD generator" group p[. 

In contrast to HERWIG and Pythia, which use backward-evolution algorithms, the GGPSl 
and GGPS2 programs evolve forwards, i.e. upwards in scale toward the hard process. This 
has the advantage that the evolution of the structure function is generated by the Monte 
Carlo algorithm itself, allowing a non-trivial test, whereas backward evolution is guaranteed 
to reproduce whatever distribution function is input. Like most implementations of forward 
evolution, GGPSl and GGPS2 produce weighted events, which can be inconvenient for detector 
simulation. Forward evolution can also be extremely inefficient, particularly at small x, but 
this problem is solved in GGPSl and GGPS2, as described in Section [4.3| . 

In the Dipole Cascade Model, implemented in the Ariadne program, there is no explicit 
initial state radiation. Instead, in DIS, all gluon emissions are described in terms of dipole 
radiation from the colour dipole between the struck quark and the photon remnant. The 
radiation is suppressed in the remnant direction because of its spatial extension. 



2.7 Hadronization 



The hadronization of the produced partons is done according to the Lund string fragmentation 
model |11]] except in HERWIG, where a cluster fragmentation model is used [0. Both these 
models are explained in detail in the report from the "QCD generator" group 0. 
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The difference in hadronization between 77 and e"'"e~ events is the presence of the photon 
remnant in the former. In Pythia the remnant is taken to be a single parton (optionally quark 
and an antiquark in case a gluon is taken out of the photon) which is treated like any other 
parton. In HERWIG, however, the cluster containing the remnant is treated a bit differently 
from the others as described in section [4.4[ 

The hadrons are subsequently decayed. Those produced by the hadronization process are 
generally taken to be unpolarized and thus decay according to pure phase-space, with standard 
decay tables as described in Ref. 0. On the other hand, when vector mesons are produced 
elastically they are strongly transverse polarized, so they are decayed accordingly, with the two 
pions predominantly taking similar energies. 



3 Comparisons 

In the following, we will compare the predictions for LEP 2 of different generators for different 
kinds of processes. Unless specified otherwise, comparisons are made on the generator level, 
using a beam energy of 87.5 GeV, a total integrated luminosity of 500 pb~^ and requiring at 
least an invariant mass of 2 GeV of particles in the region | cos6\ < 0.97. 



3.1 Exclusive channels and resonances 



Lepton pair production, e"'"e~ e^e~l^l~ is well described by the exact matrix elements, which 
are dominated by the multiperipheral diagram. In terms of the luminosity function, the lepton 
pair production is given by the QED structure function for 77 — > We have compared 

muon pair production by PC [|13| applying the luminosity function to the exact calculations 
of matrix elements by Diag36 [|14| and Vermaseren The mass threshold imposed on the 
pair isW > 300 MeV. The distributions of = max(— g^, — g|) and W obtained are shown in 
Fig. |. 



The cross-section of 77 — i? for a narrow resonance of spin- J is given by 0,|T6 



.(-,-, ^fi)=8.(2JH-l)F-( A ,|) ^^^,_;^,^^^ _^,^, (7) 

Here the resonance has mass m^, radiative width P^^, and energy dependent width P. For a 
wide resonance the Breit-Wigner term of Eq. ^ is multiplied by {mn/W)"' |]16|, with n=l,2 for 



the resonance decaying into two, three stable particles respectively. The dependence is given 
by the form factor function that satisfies -F^(0, 0) = 1. The VMD model predicts 

F^iQi^Ql) = Ti ^2^/ 2 \2 Tx ^2! 2 IZ^v-^l, V = p,u),(\),Jl'\l) ... 

(1 - gi (1 - glwj^ V 

(8) 
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Figure 1: The and W distributions of e^e — > e^e /I'^ii at ^/s — 
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Figure 2: The rjc cross-section versus beam energy: Low approximation, dashed; narrow-width 
approximation, solid; Respro results, points. 



In generating two-photon resonances, the mass spectrum is the product of the Breit-Wigner 
distribution and the two-photon invariant mass spectrum obtained from the luminosity function. 
The decay products are further described by the matrix element according to the spin-parity 
and helicity state that couple to the angular distributions of the final state particles. 
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Figure 3: T/ie 02(1320) simulations of PC in Sir phase space and TwoGEN with J^=2+, helicity 
X=2 and intermediate state pvr. Detector acceptance are set to 6 > 10°, 25° for calorimeter and 
central tracker respectively. 



We have compared the cross-section generated for rjc resonance as a function of beam energy 



which is shown in Fig. g. The dashed hne is the Low approximation [|T7| (i.e. the narrow- 
width approximation and the leading term of the EPA, but with the x-dependence of the 
integration neglected: / dP^/P^ — > logs/ml), the solid line is the luminosity function of Eq. |] 
using the approximation of W"^ = AwiW2 for small g^, and the points are the Respro [|1^ 
calculation using exact W1W2 without form factor. 



We have instrumented the photon flux generated by PC and TwoGEN [^] to simulate 
02(1320) mesons in 7r"^7r~7r° final states at LEP 1 with detector acceptance set to > 10°, 25° 
from the beam direction for calorimeter and central tracker respectively. Shown in Fig. ^ 
and b are the invariant mass and distributions of 02(1320) generated by PC and TwOGEN. 
Instrumented in PC is a narrow resonance of Eq. |^ decaying into phase space, while in 
TwoGEN the Breit-Wigner of a wide resonance is multiplied by a matrix element of J^=2~^ 
with helicity A=2 and intermediate state of 02 — > pvr — Tr+vr^vr'^. A second Breit-Wigner term 
for p suppresses events at low W. Within the detector acceptance good agreement is seen. 
The coupling to spin-parity and helicity state is demonstrated in Fig. ^ and d in which the vr" 
distribution in polar angle and acollinear azimuthal angle of the two charged pion tracks are 
compared with the phase space distributions of PC. 
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Figure 4: Transverse energy flow and transverse energy distribution for photon-photon colli- 
sions. The full (dotted) line presents the predictions by Pythia (Puojet). 




Figure 5: Pseudorapidity distribution of charged particles and multiplicity distribution for non- 
single diffractive photon-photon collisions. The full (dotted) line presents the predictions by 
Pythia (Phojet/ 

3.2 Minimum bias events 

In the following, some predictions for photon-photon collisions of the two minimum bias gener- 
ators Pythia and Phojet are compared at fixed photon-photon energy. The calculations have 
been done for inelastic, non-diffractive events, using the photon parton distribution functions 
SaS ID (Pythia) and GRV LO flOl (Phojet). The transverse energy is one of the typical 
minimum bias quantities which can be measured without collecting very high statistics. In 
Fig. ^ the transverse energy flow as function of the pseudorapidity is shown. The transverse 
energy distribution is shown in Fig. Both generators give similar predictions. 

In Fig. I^a the charged particle distribution as function of the particle pseudorapidity is 
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Figure 6: Transverse momentum distributions of charged particles predicted by Pythia (full 
line) and Phojet (dotted line). 

shown. Pythia predicts a higher charged particle density than Phojet. This can be compared 
to the experimental value dNch/drj = 1.4 ±0.1 for inelastic proton-proton collisions pO[ at 
= 24 GeV. The shapes of the charged particle multiplicity distribution calculated with 
the generators (shown in Fig. ^d) are similar, however Phojet predicts a lower value for the 
average multiplicity. 

In general Pythia gives a larger difference between 77 and pp (and 7p) collisions than 
Phojet but this discrepancy is expected to diminish once both programs have been properly 
tuned to HERA data. 

Finally, the transverse momentum distribution of charged particles is shown in Fig. |^. 
Whereas at low energies the generator results agree, Phojet predicts at ^/s = 50 GeV more 
charged particles at high transverse momentum than Pythia. 

3.3 Deep inelastic scattering 

When measuring the photon structure function in deep inelastic e7 scattering, it is of course 
necessary to be able to accurately measure the and the Bjorken-x in each event. The 
is easily obtained from the energy and angle of the scattered electron, but since the energy of 
the target photon is unknown, x must be determined from the hadronic final state. Thus one 
measures the distribution of visible hadronic mass, PVvis, and hence x^is = Q"^/ (W^vis + This 
is then converted into the distribution of x using an unfolding procedure, typically based on 



It is clear that the more of the solid angle over which we detect hadrons, the better correlated 
W^is will be with the true hadronic mass, PVtrue, so the less work the unfolding procedure will 
have to do. For now we simply define PVvis to be the total invariant mass of all the hadrons 
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Figure 7: (a) The transverse energy flow in deep inelastic 67 scattering at LEP 2 as predicted 
by different generators for x < 0.01. (h) The mean value of VFvis o^s a function of Wtme o^t 
LEP 2 as predicted by different generators. 



within the angle covered by the tracking system of a typical LEP detector, | cos6\ < 0.97, as 
done in most previous analyses, and return to this point later. 

The unfolding procedure relies heavily on the event generator's ability to correctly model 
the hadronic final state. In previous analyses (see for example Refs. p3| , p^ ), the consistency of 
the generators used for the unfolding has been checked using inclusive distributions of the data, 
such as the distributions in Q^, W^is and particle multiplicity and momentum distributions. 
The problem is that such distributions can be fitted with basically any generator by adjusting 
the number of events in each x and bin, i.e. by modifying the input parton distribution. To 
really check the generator's description of the final states, one needs to investigate less inclusive 
distributions. 

Experience from HERA shows that at small x, the distribution of transverse energy, E±, in 
the proton direction is very difficult to describe. Indeed, programs with conventional "DGLAP- 
like" initial-state parton showers cannot explain the measured E±_ perturbatively and are dom- 
inated by the non-perturbative components. On the other hand, the dipole cascade in the 
Ariadne program describes the E± flow rather well at the perturbative level, with only small 
hadronization corrections. While more detailed measurements might hope to distinguish them, 
the most recent versions of all models are in good agreement with current data. 

In Fig. I^a we show the E± flow for small-x events at LEP 2 predicted by some event 
generators. The reach in x is not as large as at HERA, but the differences are still large. In the 
forward and central regions, the differences are mostly due to the fact that HERWIG corrects 
the hardest emission to reproduce the full 0{a'^as) and 0{a^) matrix elements, while this is 
only done partly in Ariadne and not at all in Pythia. In the backward (photon remnant) 
direction, Pythia is expected to be lower than Ariadne because of the differences in the 
parton showers as explained above. HERWIG is here higher than Pythia due to the special 
treatment of the remnant fragmentation. 
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Figure 8: The angular distribution of energy flow as a function of x at LEP 1 as predicted by 
HERWIG for events with A < Q'^ < ?,Q GeV^ . The axes are such that the tagged electron is 
always on the forward, 0° side, while the photon remnant goes in the backward direction. 

It is clear that the relationship between Wvis and W^true is closely related to the energy flow, 
and in Fig. we see that the models that give higher E_[_ flow also give a stronger correlation 
between W^^s and W^true as expected. Also, all models predict a very weak correlation at large 
W^true, because more and more of the photon remnant falls outside the assumed acceptance as 
X gets smaller. This is demonstrated more clearly in Fig. ^, where the angular distribution 
of energy flow as a function of x is shown for LEP 1. The situation is even more extreme at 
LEP 2. The amount of energy falling into the central regions of the detector hardly increases 
with W, while almost all the increase is concentrated in the far forward and backward regions. 

Our assumed detector acceptance covers all but the last flve bins at each end, but we see 
that this is where most of the energy increase lies. However, all the LEP experiments have 
detectors in this region, which are used to tag electrons, and cover all but the very last bin 
of this plot. Thus if these could be used, even just to sample some of the hadronic energy in 
this region, a greatly improved W measurement could be made. Preliminary studies indicate 
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that in addition to the neutral pions, which can be measured because they decay to photons, 
a significant fraction of charged pions deposit some of their energy on the way through the 
detector. As a generator-level prescription to get some measure of how big an improvement 
this will make, we assume that neutral pions are perfectly measured, while no other hadrons 
are measured at all. In addition, we multiply all vr" energies by a factor of three to arrive at an 
estimate of the total hadronic energy in the forward region. 

In addition to extending our angular coverage, we can use additional kinematic variables 
of the final state to improve our W measurement. Defining the light-cone components of the 
particle momenta p± = E ± can write the invariant mass of the hadronic system as 



w^'= Ep^+ Ep^- - EpI 



(9) 



where i runs over all hadrons. Using energy-momentum conservation (and neglecting the vir- 
tuality of the target photon) we can get some of the terms in Eq. ^ from the scattered electron, 
to define 



Wfec = [Pe+ - Pe'+ )\Y.Pi-) ~ 



Pe 



(10) 



This is equivalent to the well-known Jacquet-Blondel idea in photoproduction pSf except for 



the inclusion of the transverse momentum component. The sum over i now runs over hadrons 
in the central region, | cos^| < 0.97, and vr^s in the forward region, 0.97 < | cos6\ < 0.9996, and 
we include an extra factor of three for the forward vr^s. 

Using Eq. |TD| has two advantages over the naive method of just using the invariant mass 
of detected hadrons. Firstly, the detector resolution enters as the product of hadronic and 
leptonic resolutions, rather than as hadronic-squared which, since the leptonic resolution is 
usually much better than the hadronic, gives improved overall resolution. Secondly, the effect 
of missing particles in the current direction ^ ~ is minimized, because they give a negligible 
contribution to E — p^. This, in addition to the inclusion of forward 7r°s, means that W^ec is 
much better correlated with the true W, as seen in Fig. Also the differences between the 
models is much smaller. 

In anomalous events, the photon remnant typically has larger transverse momentum, which 
means that in such events the correlation between W^is and PVtrue is higher, as is seen in Fig. ^d. 
The definition of anomalous is, however, not the same in all generators, and in HERWIG, where 
anomalous events are defined by having a large transverse momentum remnant, the differences 
are larger than in Pythia, where the remnant has larger transverse momentum only on average. 

Finally in double-tag events, we have more accurate information about the target photon 
energy, and hence the W^^, as is seen in Fig. |TU|. It may be possible to use such events to 



calibrate the unfolding procedure for single-tag events. However, care must be taken, as in 
double-tag events, the target photon is more off-shell and the fraction of anomalous events are 
expected to be higher. See also sections 2 and 5 of the report from the "77 Physics" [|l| working 
group for more discussions of double-tag events. 
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Figure 9: (a) The mean value ofWrec as a function ofWtrue o-i LEP 2 as predicted by different 
generators, (b) The mean value ofW^is for anomalous and normal events as a function o/VFtrue 
at LEP 2 as predicted by different generators. 
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Figure 10: Scatterplot o/PVtme i^s. Wtag O't LEP 2 as predicted by HERWIG for double-tag events 
after detector simulation. 
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Figure 11: The inclusive jet cross-section as a function of pt according to HERWIG (solid), 
Pythia (dashed) and Phojet (dotted), when all their underlying event models are turned 
off (o-), o-nd the relative changes when they are turned on (h). Errors shown come purely from 
the statistics of the Monte Carlo samples. 

3.4 High-p_L 

Next-to-leading order QCD predictions are now available for both hadron and jet distributions 
in 77 collisions. Comparisons of data with these predictions are hoped to provide additional 
constraints on the parton content of the photon, particularly the gluon. However they apply to 
the partonic final state, consisting of at most three quarks or gluons, and cannot be compared 
to data without incorporating hadronization effects. For this to be done meaningfully, it is 
important to use general purpose QCD event generators, so that the parton level is as accurate a 
representation of the calculation as possible, while the hadronization is well-constrained by other 
reactions. Indeed the problems faced in 77 are almost identical to those in 7p photoproduction 
at HERA, and we can expect that they will have been explored in detail by the time LEP2 
analysis starts. 

The hadronization effects can be broadly split into two groups: hadronization itself, and 
'underlying event' effects, the latter coming principally from twice-resolved events in which the 
two photon remnants interact with each other in addition to the main scattering. Of course this 
separation is model-dependent, as hadrons described as coming from initial-state radiation in 
one model could be described as an underlying event in another, but it is a useful guide to where 
we expect the models to be more or less reliable. Hadronization itself is expected to be largely 
process-independent so models tuned to e"'"e~ annihilation should give a reasonable description 
of data. Underlying event effects are much more poorly understood and the models are almost 
completely unconstrained at present. Their main effect is to spray additional transverse energy 
around the event, which can severely distort jet measurements, principally by adding extra 
transverse momentum to the jets. Because the jet spectrum is so rapidly falling, this increases 
the jet cross-section considerably at any given jet transverse momentum. 
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In Fig. |TT], we show the inclusive jet cross- sect ion, as a function of the jet transverse mo- 
mentum. Jets are reconstructed using the CDF cone algorithm with a cone size of 1.0 and a 
minimum transverse momentum of 3 GeV, using all hadrons within the angular acceptance, 
|cos^| < 0.97. We see reasonable agreement between Pythia, Phojet and HERWIG when 
their underlying event and multiple scattering models are turned off. However, while all of them 
predict a significant increase with multiple scattering and the underlying event, there is little 
agreement about its size. It is worth pointing out that HERWIG's soft underlying event model 
produces far too big an effect to fit HERA photoproduction data, one of the motivations for 
incorporating the JiMMY multiple hard interaction model into HERWIG instead. Amongst the 
three other models the relative effect of multiple interactions is comparable to the differences 
between the models. It is clear that these effects must be understood before an accurate jet 
measurement can be made below about 10 GeV. 

Of course the data itself can be used to study the effects and constrain the models. The jet 
energy profile is particularly sensitive to the underlying event, since perturbative radiation and 
hadronization are concentrated at the core of jets, while the underlying background is much 
more diffuse. Furthermore, since one expects the underlying event to mainly be important in 
the twice-resolved process, if we make a physical separation of direct and resolved photons we 
can test this picture. At HERA, a cut is made on x^, the fraction of the photon's light-cone 



momentum carried by the reconstructed jets ||2^. In fact the HERA experiments define this 



in dijet events from the two hardest jets, while we propose a slightly different approach for 77 
collisions: to use all the reconstructed jets, regardless of how many there are. Clearly when 
this fraction is close to 1 there can be no photon remnant, whereas when it is much smaller 
than 1 there must be one. The cut is typically set at around 0.7. For 77 collisions, we define a 
pair of momentum fractions by analogy, as 

Sparticles(-^ i Pz) 

where, following the discussion in section |3.3| , we include three times the forward 7r° energy in 
the sum over particles, which significantly improves the measurement of the denominator of 
Eq. |ll]. We define the axes such that > x~ . 

In Fig. 0, we show the Ex profile of the hardest jet in the central region, |?7| < 1, where 1] is 
the jet rapidity. We only include the transverse energy within the azimuthal region | A0| < 1, 
so in a two-jet event the Et of the other jet should not contribute to the pedestal of this jet. We 
see that for the direct events, in which both light-cone momentum fractions are large so there 
is no significant underlying event, the multiple scattering makes very little difference. On the 
other hand, for twice-resolved events, in which both photons have remnants, it makes a lot more 
difference. The jet pedestal is sufficiently raised to increase a jet's transverse momentum by 
about 200 MeV. While this may not seem a large shift in absolute terms, the jet cross-section 
actually decreases by about 25% in going from 3 to 3.2 GeV, giving rise to the correction 



predicted by Phojet in Fig. |ll|b. Because of the strong correlation between the shape of the 
jet and the shift in the cross-section, it seems hopeful that the models can be constrained by 
data, to provide a reliable unfolding of these effects. Indeed this is already in progress at HERA. 
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Figure 12: The jet profile of central jets divided into (a) direct events and (h) double-resolved 
events according to a physical definition, as predicted by Phojet with (solid) and without 
(dotted) multiple interactions. At] is the rapidity relative to the jet's, with the axes defined such 
that > x~ such that there is always more of a remnant to the left than the right. 

Of course, one would like to make a more direct measurement of the nature of the un- 
derlying event, to differentiate between soft and multiple-hard interaction models. However 
it is extremely difficult to find event features that are unambiguous in this respect, as one's 
naive picture of four jets in back-to-back pairs does not survive hadronization and realistic jet 
definition, so one is forced to look in the hard tails of distributions that are usually not well 
predicted anyway. Nevertheless, if direct evidence of multiple hard interactions could be found 
it would be extremely important for our understanding of photon and hadron collisions as a 
whole, and it is certainly worth continuing to search for such signals. 



3.5 Heavy quarks 

In this study we have compared five different generators of charm production in two-photon 



events: Vermaseren; Pythia; HERWIG; GGHVOl |23,|8[; and MiNIJET [||]. All of these 
can generate the direct process, but only the last 4 the resolved process. For the comparison 
we chose the same parameters for each model: charm mass=1.7 GeV; minimum cc invariant 
mass=4.0 GeV; beam energy=85.0 GeV; and the GRV parton distribution set (for the resolved 
process) . 

We generated 10000 events with each program and compared the following distributions: 
scattered electron energy; final state invariant mass; energy, pt and rapidity of the charm quark; 

of charm quarks with cos(^) < 0.9 and energy > 2.0 GeV. In most cases the differences 
between the generators turn out to be rather small, so only a small selection are shown here 
(Figs. p!3HT5|). Since they are so similar, we emphasize the differences by plotting the absolute 



distribution for only one generator (Vermaseren/PYTHIA for direct /resolved production) and 
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Figure 13: Comparison of five generators of the direct production of charm quarks in two 
photon collisions. The top left histogram is the distribution for the Vermaseren generator. 
Other generators are shown as a ratio to the Vermaseren distribution in remaining plots. 



showing the results for the other generators as a ratio to it. 

It can be seen that in the case of direct production the programs produce very consistent 
results. It should be noted that the Vermaseren generator is a full matrix element calculation 
while all the other generators involve a convolution of some version of the Equivalent Photon 
Approximation luminosity function with a cross-section for real photons. Thus it should be 
considered the correct answer to which the other calculations should be compared. 

For the resolved process there is a larger variation between the generators. Unfortunately at 
this stage it is not clear which if any of these should be regarded as standard. The distribution 
that most clearly shows the differences is that of the rapidity of the charm quark (Fig. |T5|), 



where the Minuet generator produces quarks that are more forward peaked, while GGHVOl 
produces more quarks in the central region. Pythia and HERWIG are fairly similar and lie 
somewhere in between. 

Experimental charm measurement is generally restricted to the central region, with rapidity 
less than about 2. In Table ^ we show the selection efficiency of an imaginary experiment with 
100% tagging efficiency for all charm quarks within the acceptance cuts, cos(6') < 0.9 and 
p_L > 2.0 GeV. The differences in the rapidity distribution for resolved events are directly 
reflected in the selection efficiencies. One can see that an experiment using Minuet to correct 
their data would 'measure' a resolved cross-section almost a factor of two larger than one 
using GGHVOl, with Pythia and HERWIG lying between the two. It is therefore clearly a 
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Figure 14: Comparison of four generators of the production of charm quarks via the single re- 
solved process in two photon collisions. The top left histogram is the distribution for the P ythia 
generator. Other generators are shown as a ratio to the Pythia distribution in remaining plots. 
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Figure 15: Comparison of four generators of the production of charm quarks via the single re- 
solved process in two photon collisions. The top left histogram is the distribution for the Pythia 
generator. Other generators are shown as a ratio to the Pythia distribution in remaining plots. 
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Generator 


Direct 


Resolved 


Vermaseren 


42 




GGHVOl 


39 


28 


Pythia 


39 


24 


HERWIG 


40 


21 


Minuet 


37 


16 



Table 2: Percentage of events passing 'acceptance' cuts. 



high priority to understand where these differences arise and in particular whether the spread 
represents a genuine uncertainty in the measurement or is simply an indication that we should 
not trust one or some of the models. Unfortunately little progress has been made with this 
during the workshop. 

Comparisons with data for the event features, rather than just the total cross-section, should 
also help to resolve this discrepancy. Within the limited statistics (33 events) Aleph found that 
GGHVOl gave a better description of data than the Vermaseren generator alone . 



4 Description of programs 

In this section we detail some of the event generators for 77 physics, concentrating mainly on 
those in which there has been significant development during this workshop. Other programs 
are commented on and described briefly in section ^. 



4.1 Ariadne 



Version: 
Author: 



Program size: 
Program location: 



4.07 of 15 August 1995 [|0| 
Leif Lonnblad 

NORDITA, Blegdamsvej 17, 
DK 2100 Copenhagen 0, Denmark 
Phone: + 45 - 35325285 
E-mail: leif@nordita.dk 
12839 lines 

[http:/ /suryall.cern.ch/users/lonnblad/ariadne/ 



The Ariadne program implements the Dipole Cascade Model (DCM) for QCD cascades 



31],^. In this model the emission of a gluon gi from a qq pair created in an e"'"e~ annihilation 
event can be described as radiation from the colour dipole between the q and q. A subsequent 
emission of a softer gluon g2 can be described as radiation from two independent colour dipoles. 
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one between the q and gi and one between gi and q. Further gluon emissions are given by three 
independent dipoles etc. In this way, the end result is a chain of dipoles, where one dipole 
connects two partons, and a gluon connects two dipoles. This is in close correspondence with 
the Lund string picture, where gluons act as kinks on a string-like field stretched between the 
qq pair. 

Further details of how Ariadne generates the QCD cascade in e"'"e~-annihilation can be 



found in Ref. |30] and elsewhere in these proceedings M. Here only the parts relevant for e7 



DIS and high-p^ 77 is presented. 

The treatment of DIS is very similar to e"*"e~ ||33| gluon emission is described in terms of 
radiation from the colour dipole stretched between the quark, struck by the electroweak probe, 
and the photon remnant. The difference is that, while q and q are both point-like in the case of 
e"'"e~, the photon remnant in DIS is an extended object. This results in an extra suppression 
of radiation in the remnant direction. 

The suppression depends on the transverse size of the remnant, which is taken to be inversely 
proportional to its intrinsic transverse momentum ki±. In anomalous events where ki± is larger, 
the suppression is therefore smaller than in VMD-like events. Also in events where the target 
photon is significantly off-shell, the transverse size is taken to be the inverse of the maximum of 
ki± and the photon virtuality. Similarly in the case of events with <C W"^, where gluons may 
be radiated with transverse momentum larger than Q, also the struck quark may be treated as 
extended with a transverse size oc 1/Q. 

The DCM is evidently very different from conventional initial-state parton shower models. 
Tracing emissions "backwards" from the struck quark as in an initial-state shower, these would 
be unordered in transverse momentum, while in a program like Pythia these emissions would 
be ordered in falling transverse momentum. In the DCM, the transverse momentum of gluons 
far away from the struck quark, i.e. close to the remnant, can be much larger than in an 
initial-state parton shower. 

The first emission in the dipole cascade is corrected to match the 0{as) matrix element for 
gluon emission. The gamma-gluon fusion diagram is, however, not yet included for the 67 DIS, 
because of technical difficulties in the current interface to Pythia. 

The interface to Pythia, described in the report from the QCD generator group also 
enables Ariadne to generate high-p^ 77 scattering. Here, Pythia is used to generate the 
hard sub-process. The outgoing partons and remnants are then connected with dipoles, which 
are allowed to radiate, restricting the transverse momentum of the emissions to be smaller than 
that of the hard interaction and treating all remnants as extended objects as in DIS. 

For 77 it is also possible to run Ariadne with Pythia using multiple interactions. However, 
this part of the interface is still preliminary, and more studies are needed. 
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4.2 GGHVOl 



Version: Version 1.0 of 1/11/1994 

Author: M. Kramer, P. Zerwas (DESY), 

J. Zunft (formerly of DESY, no longer in HEP) 

A. Finch, 

School of Physics and Materials, 
University of Lancaster, 
Lancaster LAI 4YB United Kingdom 
Phone: (+44) 1524 593618 
E-mail: A.Finch@lancaster.ac.uk 

Program size: Interface routines: 1148 lines 

Physics routines: 2639 lines 
Integration package (BASES): 4432 lines 

Program location: ftp://lavhep.lancs.ac.uk/ 



and 



GGHVOl was developed by the above authors as a Monte Carlo implementation of the 
calculation in Ref. [^] of heavy flavour quark production in gamma gamma collisions at next 
to leading order. It can either produce direct or single resolved events at a range of beam 
energies, but is restricted to real photons only. 

The program is not a complete NLO generator. It does generate the 2^3 process according 
to the NLO matrix elements, which is cut off against soft and coUinear divergences. The 
remaining events are, however, approximated with 2^2 processes generated only to leading 
order, but rescaled so that the total NLO cross-section is reproduced. 

During the initialization stage the routine FANDK calculates the factor by which the 2^2 
process should be increased to achieve the correct cross-section. In order to do this it needs 
to know beforehand the total cross-section. This it finds from a parameterization of the total 



cross-section from Ref. ||2^. This is a quick method but less flexible than recalculating the 
total cross-section from scratch. The latter approach may be adopted in a later version. The 
chief criticism of this approach is that there is just one global correction factor whereas a more 
sophisticated approach would be to calculate a different factor for different regions of phase 
space. 



The main event generation loop uses the coupled routines BASES and SPRING |0 which 
together constitute a general purpose integration and event generation package. The input to 
both routines is a function to be integrated as a function of n random variables. In this case the 
routine GENHVY. This routine uses the input random numbers to first pick whether to generate 
a 2 or 3 body event, it then calculates the fraction of the energy of the incoming beam electrons 
taken by each photon. In the resolved photon case the fraction of the photon's energy taken by 
a gluon is also found using the Gliick Reya and Vogt next to leading order parton distribution 



functions [10|. Momenta for the outgoing particles are then chosen using the subroutine RAMBO. 



The weight for the event is found using routines based on the NLO calculation of Ref. 
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multiplied by the phase space weight from RAMBO. The 2 body ('Born term') contributions 
is additionally scaled by the factor calculated in the initialization step as described above. 

During the BASES stage the program first optimizes a 'map' of the function to provide effi- 
cient calculation, and then calculates the total cross-section (for GGHVOl this simply provides 
a check as the total cross-section is already known). During the SPRING stage the function 
is repeatedly called to generate unit weight events which are put in the LUJETS common for 
fragmentation by the LUEXEC routine from the Jetset package. The map generated by BASES 
is written out to a file, which means the program can later be run from the SPRING step only. 
In this case it reads back the map produced in a previous run and by simply varying the random 
number seed a fast event generation can be achieved.. 



4.3 GGPSl, GGPS2 



Version: 1.0 05 



Author: T. Munehisa, K. Kato, D. Perret-Gallix 

Phone: +81-552208584 (TM), +81-333421264 (KK), 
+33-50091600, +41 22 7676293 (DPG) 
E-mail: munehisa@hep . esb . yamanashi. ac . j p , 
kato@sin.kogakuin.ac.jp, perretg@cernvm.cern.ch 
Program size: 2213 (GGPSl) and 2841 (GGPS2) lines 
Program location: ftp:/ /lapphp8. in2p3.fr/pub/keklapp/ggps/ggps.tar.gz 

This generator simulates jet productions in two-photon process based on the leading-log (LL) 
parton shower (PS) technique [^]. Two cases are separately treated, namely, the deep inelastic 
scattering of the photon (GGPSl) and the scattering of two quasi- real photons (GGPS2) ||35|| . 
Both processes begin by the PS space-like evolution, then the hard scattering of partons takes 
place, followed by the time-like evolution of the final state partons. 



The non-singlet quark distribution in the photon, qNs{x,Q^), obeys 

(12) 

The inhomogeneous term, k^^si^)^ proportional to + (1 — x)^. a is the QED coupling con- 
stant, the QCD coupling constant is defined as: as{Q'^) = 47r/(/?o InQ^/A^) = ao / (In Q"^ / A'^) 
with /?o = 11 — {2/3)Nf- Nf is the number of fiavours. Eq. |12| can be brought into the integral 
equation 



qNs{x,Q')= f'^-^K^^l{x/y,s)qr,s{y.Ql)+rl%f^-^^^^^ (13) 
Jx y JQl K"^ Jx y iTi 



where s = ln(a,(Q2)/a^(Q2)) ^ ln(ln(QVA2)/ ln(Q2/A2)) and r] = ln(ln(QVA2)/ In(irVA')). 
The first term represents the vector meson dominant part (VMD) while the second corresponds 
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to the perturbative photon. Here K^^g{x,s) is the QCD kernel function defined by the 
inverse MeUin transformation 

A-<»^(...)=r^-*.-»e"*)', (14) 

JrQ—ioo ATTl 

where d{n) is the moment of Pj^^j^x) / {2tt) . 

The singlet distribution is handled in the similar manner. 



Due to the inhomogeneous term in Eq. |12|, the total energy is conserved only if the photon 
energy is included, 



rfa:xE(g/(x,Q^)+g;(x,Q2)) + G(a:,Q^)] = 
In the generation the right-hand side is used as the weight of the event. 



Let us summarize the main steps present in the algorithm generating the partons in the 
course of the evolution. 

1. Selection of a Q^. 

2. Calculation of the energy of the VMD part (independent of Q"^) and of the energy of the 



perturbative photon part by using Eq. [15|. The sum of these energies is used as the weight 
of the event. 

3. Selection of the actual process, either VMD or perturbative photon, according to the ratio 
of energies. 

4. If VMD is chosen, the usual QCD evolution from (the cut-off momentum) to is 
performed. 

5. In the case of a perturbative photon, the virtual mass squared according to the 
probability dK'^ / K'^ is determined. Then the flavour of the partons are selected according 
to the ratio of charges squared. 

6. For each quark or anti-quark, the usual QCD evolution from up to is performed. 

This algorithm is common to the deep inelastic photon scattering and to the two quasi-real 
photon scattering with large The hard scattering parts are, however, different. In quasi-real 
photon collision the initial states g-g, g-g, g-G, G-G, g-7 (7-7 and 7-7 are taken into account for 
the hard scattering, although in the virtual-quasi-real case only the 7*-g and 7*-G' subprocesses 
are considered. 

The reference cross-section used to select events is the differential cross- sect ion: 

da^/dpl = n/p^. (16) 
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The ratio of the hard cross-section to Eq. |T6| is counted as the weight of the event. The 
argument of in the hard cross-section is set to be p^, because there are some processes 
whose cross-section is dominated by the t-channel contribution. After the time-hke evolution 
of the produced partons has been performed, the energy of the initial photons is fixed and thus 
the four-momenta of all partons are determined. 

Events are generated with a weight whose maximum is unknown at the beginning of a run. 
A maximum value is arbitrary fixed by the user, overweight events are rejected and counted. At 
the end of the run, if the number of overweight events is too large, compromising the generation 
accuracy, the user must increase the maximum weight, reducing consequently the generation 
efficiency. 



Finally the generated partons are hadronized following the Jetset |^ mechanism. Colour 
flows are properly matched and adopt the Jetset definition. Two types of colour singlet exist: 
a string beginning and ending with quarks possibly embedding one or many gluons or a closed 
colour loop of gluons. 



4.4 HERWIG 



Version: 
Authors: 



Program size: 



5.8d of October 1995 || (and 5.9 of December 1995) 

G. Marchesini^ B.R. Webber^, G. Abbiendi^, I.G. Knowles^ 

M.H. Seymour^, L. Stanco^ 

^Dipartimento di Fisica, Universita di Milano. 

^Cavendish Laboratory, University of Cambridge. 

^Dipartmento di Fisica, Universita di Padova. 

^Department of Physics and Astronomy, University of Glasgow. 

^Theory Division, CERN. 

E-mail: webber@hep.phy.cam.ac.uk, knowles@v6.ph.gla.ac.uk, 
sey mour @sur ya 1 1 . cern . ch . 
15500 lines 



Program location: http://suryall.cern.ch/users/seymour/herwig/ 



HERWIG is a general-purpose QCD Monte Carlo event generator for simulating Hadron 
Emission Reactions, With Interfering Gluons. Its general design philosophy is to provide as 
complete as possible an implementation of perturbative QCD, combined with as simple as possi- 
ble a model of non-perturbative QCD. It does this uniformly for a very wide range of processes, 
allowing the parameters to be fitted in one reaction, principally e"''e~ annihilation, and applied 
to other reactions. Although HERWIG has been capable of simulating 77 collisions for some 
time, there were many technical deficiencies, which practically prevented experiments from us- 
ing it. Many of these have been rectified during the workshop, resulting in the preliminary 
versions 5.8a, 5.8c and 5.8d. However work is still ongoing, and some of the features described 
below will not be available until the full version release towards the end of this year. Since 
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much of HERWIG is described in Ref. 0, we concentrate on the additional features relevant to 
77 collisions here. Original references can be found in Ref. ^ . 



Event generation proceeds much like the general case described in the introduction. The 
equivalent photon approximation (EPA) is used, correctly generating the dependence. Since 
this means that the photons no longer collide head-on, they are boosted to their centre-of-mass 
frame, where the remainder of the event generation is performed. At the end of the event they 
are automatically boosted back to the lab frame. 

Since the user can control the range generated (through the variables Q2WWMN and 
Q2WWMX), it is in principle possible to generate all event classes using the EPA. However, 
to provide a more accurate description of high-Q^ tagged events, these are described as deep 
inelastic scattering, e7 — > e+hadrons including the full electron kinematics, rather than as 
7*7 — i>hadrons using the EPA. 

The treatment of deep inelastic lepton-photon scattering is essentially identical to that of 
lepton-hadron, except for the inclusion of the point-like photon-quark coupling in the initial- 
state parton shower. The hard process is generated as eq eq according to whichever parton 
distribution function is selected from PDFLIB. This is controlled by the variables AUTPDF and 
MODPDF, which hold the 'author group' and set number respectively, for example GRVph and 3 
for the leading-order photon set of Gliick, Reya and Vogt. 

The outgoing quark produces a parton shower exactly like that in the final state of e^e^ 
annihilation, described in Ref. 0. The incoming quark also produces a parton shower, which 
is generated using the 'backward evolution' algorithm. One can imagine this as an evolution 
in the factorization scale from the large scale of the hard process down towards zero — as it is 
reduced, more and more radiation is resolved, i.e. transferred from the evolution of the parton 
distribution function to the coefficient function. The inhomogeneous term in the evolution 
equation corresponds directly to the addition of a 7 — >■ qq vertex, which is straightforwardly 
included by considering the photon to be an additional parton type with a delta-function 
distribution. This results in a dynamic separation of events into point-like and hadronic. At 
small and medium x values (x^O.3), this separation is similar to that made in the input 
distribution functions, demonstrating the self-consistency of the backward evolution algorithm, 
but at large x it becomes increasingly different from them, classifying almost all large-x events 
as hadronic, whereas most distribution function sets classify them as point-like. The difference 
can be traced to the fact that HERWIG uses a cutoff in transverse momentum to separate the 
perturbative and non-perturbative regions, whereas the distribution functions make a cut on 
the virtuality of the internal line. At large x, the two differ by a factor of (1 — x) giving large 
differences in the separation scales. While it could be argued that transverse momentum is a 
more physical scale to use, as in the FKP model, a more conservative approach is to simply 
say that this indicates a region of uncertainty, and any analysis that relies on the classification 
of this large x region into the two components should be considered highly model-dependent. 

All partons produced by the initial-state cascade undergoing further parton showering just 
like any other outgoing parton. 
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Colour coherence is as important in initial-state cascades as it is in the final state. However 
it is less well understood and only an approximate treatment is implemented in HERWIG. At 
large x, all emitted gluons must be soft, the emission pattern is identical to a final state dipole, 
and the emission is angular-ordered. At small x, rather than restricting soft gluon momenta, 
we actually 'look inside' the soft gluons, resulting in a different coherence structure. At 'fairly 
small' X, this can be approximated by imposing transverse momentum ordering, which is how it 
is implemented in HERWIG, but at very small x the correct approach is to use angular ordering 
and modified emission probabilities. Although this has never been implemented as a full Monte 



Carlo event generator, a partial implementation was discussed in Ref. and shown to give 
similar results to the HERWIG algorithm at the x values expected at LEP 2 or even HERA. 
The actual evolution variable used by HERWIG smoothly interpolates the two regions. 

Just as in e+e^ annihilation, HERWIG is not capable of covering the whole of phase-space 
with its parton shower emission. This is corrected using the methods proposed in Ref. | 4^ , 



which ensure that the hardest emission (which is not necessarily the first owing to the ordering 
of opening angles) agrees with the exact matrix element (the sum of higher-order resolved, 
g qq, q qg, and point-like, 7 qq). It is worth noting that azimuthal correlations 0] 



are correctly included within the 'dead- zone' region of Xp and Zp defined in Ref. |^3|, but not 
within the parton shower itself. 

After the parton cascade, the system is hadronized according to the cluster model. For 
outgoing partons, this is identical to e"'"e~ annihilation, described in Ref. 0, but for hadronic 
events, the photon remnant is treated specially. It is given a limited transverse momentum of 
width PTRMS, and the cluster that contains it is given a Gaussian mass-squared distribution, 
resulting in a limited pt spectrum of produced hadrons, as expected for such a soft object. An 
'underlying event' model is provided, which increases the energy released during the break- 
up of the remnant, but it seriously overestimates HERA data on DIS at small x, and it is 
recommended that it be switched off by running process IPR0C=19000 rather than 90000. 

Untagged and low-Q^ tagged events are generated using the equivalent photon approxima- 
tion to split both beams to photons. One hard process type is selected from those hsted in 
Ref. [0. At present there is no facility to mix events of different types in the same run, although 
it is clear that this would improve the utility of the program and it is a planned improvement. 
Events are generated according to the leading order matrix elements for 2 — >■ 2 scattering, and 
initial-state and final-state showers are added exactly as in deep inelastic scattering. Another 
current problem is that the phase-space cuts accessible to the user are completely different for 
direct and resolved processes, making it difficult to generate both uniformly. Fixing this is 
another planned improvement. 

In all event classes except deep inelastic scattering, it is recommended that the soft un- 
derlying event be selected. This models the collision between the photon remnants as a soft 
hadron-hadron scattering, which produces a uniform rapidity plateau of extra hadrons on top 



^This applies to versions 5.8d onwards. In earlier versions, the two processes would be expected to bracket 
the correct result. 
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of those from the perturbative event. This is based on the UA5 minimum bias model, and is 
essentially just a parametrization of data. The soft VMD scattering process, IPR0C=8000, uses 
the Donnachie-Landshoff cross-section and generates events as in the soft underlying event. 

As discussed earlier, one would expect multiple hard scattering to be an important effect 
in untagged 77 collisions at LEP 2. Its contribution to the underlying event is also a major 
source of uncertainty in the measurement of jets at HERA and it is important to have several 
different models of it, to get some estimate of the uncertainty in the predictions. Although 
such scattering is not included in HERWIG at present, there is an interface to the Jimmy 
Generator (which can also be obtained from the web-page listed above). However, there is 
not yet a smooth transition between hard and soft multiple scattering, and Jimmy and s.u.e. 
should be considered mutually exclusive at present. 

The default parameter set that comes with HERWIG is tuned to OPAL data in the case of 
parameters that affect e"'"e~ annihilation, but are simply theoretically prejudiced guesses for the 
others. At present the HERA experiments are involved with tuning the additional parameters 
and these will be included in future releases, hopefully improving the predictivity for 77 physics. 



4.5 PHOJET 



Version: 
Author: 



Program size: 
Program location: 



1.04 of 20 October 1995 
Ralph Engel 

Institute of Theoretical Physics 
University Leipzig 

Augustusplatz 10, D-04109 Leipzig, Germany 
Phone: + 49 - 341 - 97 32444 
E-mail: eng@tph200.physik.uni-leipzig.de 
31000 lines 

http: // www.physik.uni-leipzig. de / "engel / phojet .html 



Phojet is a minimum bias event generator for hadronic pp, 7p and 77 interactions. The 
interactions are described within the Dual Parton Model (DPM) |^ in terms of reggeon and 
pomeron exchanges. The realization of the DPM with a hard and a soft component in Phojet 



is similar to the event generator Dtujet-93 [47, 48|. Regge arguments are combined with 



perturbative QCD to get an almost complete description of the leading event characteristics. 
Special emphasis is taken on diffractive and soft interactions. Soft and hard interactions are 
unitarized together leading to the possibility to have multiple soft and hard scatterings in one 
event. 



In the following, some comments on LEP 2 specific aspects are given. In the model the 
dual nature of the photon is taken into account by considering the physical photon state as a 
superposition of a "bare photon" and virtual hadronic states having the same quantum numbers 
as the photon. Two generic hadronic states \qq) and \qq*) have been introduced to describe the 
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hadronic piece of the photon. The low-mass state \qq) corresponds to the superposition of the 
vector mesons p, uj and (j) and a tt+tt^ background. The state \qq*) is used as an approximation 
for hadronic states with higher masses. The physical photon reads 

|7) = |7bare) + ^ + ^ W) ■ (17) 

Jqq Jqq* 

The interaction of the hadronic states via pomeron/reggeon exchange is subdivided into pro- 
cesses involving only soft processes and all the other processes with at least one large momentum 
transfer [hard processes) by applying a transverse momentum cutoff p™^°^ to the partons. On 
Born-graph level, for example, the photon-photon cross-sections is built up by: (i) soft reggeon 
and pomeron exchange, (ii) hard resolved photon-photon interaction, (iii) single direct inter- 
actions, and (iv) double direct interactions. The soft pomeron cross-sections is parametrized 
using Regge theory. The hard cross-sections are calculated within the QCD Parton Model 
using lowest order matrix elements. For soft processes, photon-hadron duality is assumed. The 
energy-dependence of the reggeon and pomeron amplitudes is assumed to be the same for all 
hadronic processes. Therefore, data on hadron-hadron and photon-hadron cross-sections can 
be used to determine the parameters necessary to describe soft photon-photon interactions. 

The amplitudes corresponding to the one-pomeron exchange between the hadronic fluc- 
tuations are unitarized applying a two-channel eikonal formalism similar to Ref. The 



probabilities / fqq and 6^//^^ to find a photon in one of the generic hadronic states, the cou- 
pling constants to the reggeon and pomeron, and the effective reggeon and pomeron intercepts 
cannot be determined by basic principles. These quantities are treated as free parameters and 
determined by cross-section fits Once the parameters are fitted, the model allows for 

predictions on photon-photon collisions without new parameters. 

The probabilities for the different partonic final state configurations are calculated from the 
discontinuity of the scattering amplitude (optical theorem). Using the Abramovski-Gribov- 
Kancheli cutting rules PU[ the cross-section for graphs with soft pomeron cuts, Ic hard 



pomeron cuts, mc triple- or loop-pomeron cuts, and ric double-pomeron are estimated. For 
pomeron cuts involving a hard scattering, the complete parton kinematics and flavours/colours 
are sampled according to the Parton Model using a method similar to Ref. extended to 
direct processes. For pomeron cuts involving parton configurations without a large momentum 
transfer, the partonic interpretation of the Dual Parton Model is used: photons or mesons are 
split into a quark-antiquark pair whereas baryons are approximated by a quark-diquark pair. 
The longitudinal momentum fractions of the soft partons are given by Regge asymptotics ||51|| . 
One obtains for the valence quark [x) and diquark (1 — a;) distribution inside the proton 
p{x) ~ (1 — xY'^ I \fx and for the quark antiquark distribution inside the photon p[x) ~ 

1 / \jx{\ — x). For multiple interaction events, the sea quark momenta are sampled from a 
p[x) ~ \jx distribution. The transverse momenta of the soft partons are sampled from an 
exponential distribution in order to get a smooth transition between the transverse momentum 
distributions of the soft constituents and the hard scattered partons. 
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In difFraction dissociation or double-pomeron scattering, the parton configurations are gen- 
erated using the same ideas described above apphed to pomeron- photon/pomeron scattering 
processes. According to the kinematics of the triple- or loop-pomeron graphs, the mass of the 
diffractively dissociating systems is sampled from a distribution. The momentum 

transfer in diffraction is obtained from an exponential distribution with mass- dependent slope 
(see Ref. [^). For the parton distributions of the pomeron, the CKMT parametrization with 
a hard gluonic component is used. The low-mass part of diffraction dissociation is approx- 



imated by the superposition of high-mass vector mesons. In order to take into account the 
transverse polarization of quasi-elastically produced vector mesons, diffractively scattered p, u 
and are decayed anisotropically. 

Finally, the fragmentation of the sampled partonic final states is done by forming colour 
neutral strings between the partons according to the colour fiow. In the limit of many colons in 
QCD, this leads to the two-chain configuration characterizing a cut pomeron and a one-chain 
system for a cut reggeon. In hard interactions the colour fiow is taken from the matrix elements 
directly [^. The leading contributions of the matrix elements give a two-chain structure which 
corresponds to a cut pomeron. The chains are fragmented using the Lund fragmentation code 
Jetset 7.3 i3l 



In order to get the LEP 2 kinematics, the complete lepton-photon vertex for transversally 
polarized photons is simulated in the program. The lepton (anti-) tagging conditions can 
be specified by the user. However, in the model only photons with very low virtualities are 
considered at the moment. The extension to virtual (and longitudinally polarized) photons is 
in progress. 

An example input file for a 77 run can be in found in the Phojet manual [5^]. 



4.6 PYTHIA/ JETSET 

Versions: Pythia 5.720 of 29 November 1995 

Jetset 7.408 of 23 August 1995 

Author: Torbjorn Sjostrand 

Department of Theoretical Physics 
University of Lund 

Solvegatan 14A, S-223 62 Lund, Sweden 
Phone: + 46 - 46 - 222 48 16 
E-mail: torbjorn@thep.lu.se 

Program size: 19936 + 11541 lines 

Program location: http://thep.lu.se/tf2/staff/torbjorn/ 



Pythia/Jetset is a general-purpose generator of high-energy particle physics reactions. 
An introduction and a survey of processes of interest for LEP 2 physics are found in the QCD 
generators section, while the full description is in Ref. [^. Here only aspects specific to 77 
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physics will be summarized. These have been developed together with Gerhard Schuler, and 
are described extensively elsewhere |]55|, p^] . 

To first approximation, the photon is a point-like particle. However, quantum mechanically, 
it may fluctuate into a (charged) fermion-antifermion pair. The fluctuations 7 ^ qq can 
interact strongly and therefore turn out to be responsible for the major part of the 77 total 
cross-section. The spectrum of fluctuations may be split into a low-virtuality and a high- 
virtuality part. The former part can be approximated by a sum over low- mass vector-meson 
states, customarily restricted to the lowest-lying vector multiplet. Phenomeno logically, this 
Vector Meson Dominance (VMD) ansatz turns out to be very successful in describing a host 
of data. The high-virtuality part, on the other hand, should be in a perturbatively calculable 
domain. Based on the above separation, Pythia distinguishes three main classes of interacting 
photons: direct, VMD and anomalous. 

For a 77 event, there are therefore three times three event classes. By symmetry, the 'off- 
diagonal' combinations appear pairwise, so the number of distinct classes is only six. These 
are: 

1. VMD X VMD: both photons turn into hadrons, and the processes are therefore the same 
as allowed in hadron-hadron collisions. 

2. VMD X direct: a bare photon interacts with the partons of the VMD photon. 

3. VMD X anomalous: the anomalous photon perturbatively branches into a qq pair, and one 
of these (or a daughter parton thereof) interacts with a parton from the VMD photon. 

4. DirectX direct: the two photons directly give a quark pair, 77 — > qq. Also lepton pair 
production is allowed, 77 i^i^ , but will not be considered here. 

5. DirectX anomalous: the anomalous photon perturbatively branches into a qq pair, and 
one of these (or a daughter parton thereof) directly interacts with the other photon. 

6. Anomalous X anomalous: both photons perturbatively branch into qq pairs, and subse- 
quently one parton from each photon undergoes a hard interaction. 

In a complete framework, there would be no sharp borders between the six above classes, 
but rather fairly smooth transition regions that interpolate between the extreme behaviours. 
However, at our current level of understanding, we do not know how to do this, and therefore 
push our ignorance into parameters of the model. From a practical point of view, the sharp 
borders on the parton level are smeared out by parton showers and hadronization. Any Monte 
Carlo event sample intended to catch a border region would actually consist of a mixture of 
the three extreme scenarios, and therefore indeed be intermediate. 

The main parton-level processes that occur in the above classes are: 

• The 'direct' processes 77 qq only occur in class 4. 

• The '1-resolved' processes 7q qg and 7g qq occur in classes 2 and 5. 

• The '2-resolved' processes qq' — * qq' (where q' may also represent an antiquark), qq — 
q'q', qq — > gg, qg qg, gg qq and gg gg occur in classes 1, 3 and 6. 

• Elastic, diffractive and low-p± events occur in class 1. 
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The notation direct, 1-resolved and 2-resolved is the conventional subdivision of 77 interactions. 
The rest is then called 'soft-VMD'. The subdivision in Pythia is an attempt to be more precise 
and internally consistent than the conventional classes allow. 

The cross-sections for the above components are based on a number of considerations. The 
VMDxVMD ones are derived from a Regge theory ansatz, with a pomeron plus reggeon form 



for the total cross-section, plus parametrizations of elastic and diffractive components ||57 |. 
The other five cross-sections are obtained by a direct integration of parton-level cross-sections 
above lower cut-offs /cq ~ 0.6 GeV for 7 qq fluctuations and pl'^"^ = 0.7 + 0.171og^(l + 
i?cm/20) [GeV] for QCD processes in the anomalous sector. The latter is a purely phenomeno- 
logical fit based on consistency arguments in the 7p sector. It does not include possible eikonal- 
ization effects, and would therefore change in a more complete treatment (studies under way). 
Taken all together, one obtains a reasonable description of the total 77 cross-section. 

The program comes with several parton-distribution sets for the photon. The default is 
SaS ID P]. In view of the above event classification, the SaS sets have the advantage that the 
separation into VMD and anomalous parts is explicit. 

Currently only real incoming photons are considered. Eventually Pythia should be ex- 
tended to (moderately) virtual photons. A separate treatment exists of the DIS region, e7 
eX, i.e. with one real and one very virtual 7, but this option is less well developed (especially 
for parton showers). Furthermore, the program does not yet generate the spectrum of real 
photons internally, i.e. it is easiest to run for a fixed energy of the 77 system. It is possible 
to use the varying-energy and weighted-events options of the program to obtain a reasonable 
photon spectrum, however. 

Some main switches and parameters of interest for 77 physics are: 

• MSEL selects allowed processes; a change to 2 would include also elastic and diffractive 
VMD events. 

• MSTP(14) sets the 77 event class among the six possibilities listed above. The most 
useful option is MSTP(14) = 10, where the six classes are automatically mixed in the 
proper proportions. Note that some variables, such as CKIN(3), differ between the event 
classes and therefore automatically are reset internally. 

• MSTI(9) tells which event class the current event belongs to. 

• MSTP(55) and MSTP(56) gives choice of parton-distribution set and library for the photon. 

• PARP(15) is the ko scale, i.e. minimum p± for 7 <-> qq fluctuations. 

• CKIN(3) sets p^min scale of hard interactions; for MSTP(14) = 10 only to be set when 
studying high-p^ jet production. 

• PARP(81) (alternatively PARP(82), depending on MSTP(82) value) sets the p^^i^ scale of 
multiple interactions in VMDxVMD events. 
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5 Other generators 



Most of the general-purpose QCD event generators described above have only become available 
for 77 physics during the last year or two. Therefore most existing analyses have used generators 
written specifically for 77 physics with little or no overlap with other processes. This means 
that they are extremely hard to test thoroughly, except on the very data they are used to 
unfold. Although of course the different programs have different strengths and weaknesses, 
they are mainly based on the same general design. In this section we describe it and make a 
few general comments on when and why it is expected to be adequate or inadequate. 

The principle way in which the older programs differ from the QCD generators is the lack 
of a parton shower stage. A hard scattering is typically generated according to the leading 
order matrix element, often taking great care with their accuracy and including many effects 
that are not yet included in the QCD programs. The final state described by these matrix 
elements ranges from a simple quark- ant iquark pair in the case of DIS to up to two jets and 
two remnants in the case of high-pf scattering. These simple partonic states are then given 
directly to the Jetset program, which implements the Lund string fragmentation model. 

It may at first sight appear that parton showers are a luxurious frill that can be added to 
the model to slightly improve event simulation, for example by including the small fraction 
of events in which a third jet accompanies the two hard jets in an event. However in the 
modern view of hadronization, they are an essential precursor to the confinement of partons 
into hadrons, which does not take place globally between the main few hard partons in an event 
but locally between nearby soft partons. It is only by adopting the latter view that one can have 
any expectation that the hadronization process is universal. Quite apart from these theoretical 
issues, as a practical issue it is certainly the case that the parameters of the string model 
need to be retuned to fit to data at different energies and in different reactions when it is not 
preceded by a parton shower, but that coupled with a parton shower a reasonable description 
of all current data can be achieved with a single parameter set. This alone is enough to mean 
that any model that uses string fragmentation without a parton shower should be considered 
descriptive rather than predictive. 

In addition to their role in setting the initial conditions for the hadronization process, parton 
showers are crucial for determining certain event features, for example the hadronic final state 
of DIS at small x, where the predictions of a 'QPM model', i.e. hard matrix element plus string 



fragmentation were found to give an extremely poor description of the HERA data ||58 



Another closely related difference between the two groups of programs is in the assumptions 
made about the photon structure function. The older programs generally make an FKP-like 
separation into the point-like and hadronic parts, a scheme with strong physical motivation, 
but then neglect the QCD evolution in the hadronic part. While one might suppose that this is 
unimportant if we only use the event generator as an unfolding tool, this is not in fact the case 
because structure functions do not just evolve by themselves, they do so by emitting gluons, 
which have an effect on the structure of the hadronic final state. QCD not only predicts that 
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changing changes the structure function, but also the features of the hadronic final state. 
This should therefore be included in any model used to unfold data. 

One possible solution would be to incorporate parton showers into the existing programs. 
However this is far from straightforward, because the evolution of a shower is strongly dependent 
on its initial conditions, namely on how the hard partons were formed. For example in DIS, 
simply setting up a qq pair and calling Jetset with its final state parton shower option switched 
on will result in emission with transverse momenta up to producing far too much hadronic 
activity. The correct solution, roughly speaking, is that the current jet produces a shower with 
upper scale of order Q, the photon remnant produces one of order its pt and the internal line 
produces an initial-state shower. There are additional complications in the high-pj case as there 
are contributions that contribute to different event classes depending on their kinematics, so 
one must impose additional constraints on the parton shower in order to avoid double-counting. 
Eventually, one realizes that to build in all the relevant conditions, the parton shower has to 
know so much about the hard process that one ends up almost having to write ones own parton 
shower algorithm from scratch. 

Should we therefore conclude that existing programs are wrong and existing analyses need to 
be redone? No, because in the kinematic regions explored so far they have several advantages. 
Firstly the range in W has been rather limited so that the move into regions in which parton 
showers and QCD effects become important has been limited. Secondly for reasonably low 
most of the hadronic event is actually seen in the detector, so the actual reliance on the models 
is rather small. Thirdly, this means that the models can be well constrained and tested by 
fitting to the majority of the event so that they are fairly predictive for the unseen remainder 
of it. Finally, and very importantly, the low energy range has many problems associated with 
it that are rather carefully treated in the dedicated generators, such as exact kinematics, target 
mass effects, higher order terms in the EPA and other polarization states of the photon. This is 
in contrast with the QCD programs, which were traditionally oriented towards the high energy 
limit and are only now starting to 'catch up' in their treatment of effects important at lower 
energy. 

In conclusion, we would say that at pre-LEP energies existing programs are perfectly ad- 
equate. At LEP 1, they have proved sufficient for most tasks, but problems are beginning 
to become apparent. For example, in DIS if one studies the energy flow in the detector as a 
function of Xvis, similar to what is shown in Fig. one obtains predictions that are reasonably 
insensitive to the shape of the structure function and hence can make a fairly stringent test 
of the model. Preliminary results seem to indicate that existing models already have prob- 



lems [^. One can fairly confidently predict that at LEP 2 these problems will be sufficiently 
severe that using the QCD models will become essential, both because of the huge leap to 
smaller x and higher and simply because of the higher statistical precision that will require 
better control of systematic errors. 

Below we give a few basic facts about some of the programs used at present. 
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5.1 DIAG36 



Author: F.A. Berends, P.H. Daverveldt and R. Kleiss |jT4 

E-mail: t30@nikhef.nl 
Program size: 4335 lines 
Program location: CPC program library 

Diag36 is an event generator for the full set of QED diagrams for e^e~ — > e"'"e~ff , including all 
fermion masses. Phase-space generation can cover the whole kinematically-allowed phase-space, 
or cuts can be applied for typical single- or double-tagged configurations. 



5.2 MINUET 



Author: A. Miyamoto and H. Hayashii p9 

E-mail: hayashii @nar as 1 . kek . j p 

The Minuet program generates the direct and the resolved photon events in two-photon 
processes of e^e~ collisions. It calculates the cross-section of light- and heavy-quark production 
according to the EPA approximation for the photon flux and the leading-order matrix elements 
for the sub-process cross-sections. It uses the BASES program for the cross-section calculation 
and SPRING for the generation of four-momenta of final state partons. The generated parton 
are hadronized using the Jetset program. It has been used in the analysis of two-photon data 
by the TOPAZ group and in background studies for the JLC Multiple interactions are not 
included in the program. 



5.3 PC 

Author: F.L. Linde iTBI 



E-mail: z66@nikhef.nl 

The PC program generates two-photon events according to Eq. ^ It provides event generation 
of fermion pairs and resonance states in phase space of up to 4 final state particles with many 
dominant resonant states. The form factor is chosen between p, vr, and J/ip poles. It provides 
a cross-section calculation of all terms in Eq. 0. 
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5.4 RESPRO 



Author: E.R. Boudinov and M.V. Shevlyagin |]T 

E-mail: BOUDINOV@vxcern.cern.ch 



The Respro program generates exphcitly the two-photon resonance in a very efficient way and 
calculates the cross-section according to Eqs. H and |^. 



5.5 TWOGAM 

Authors: S. Nova, A. Olshevski, T. Todorov [|] 

E-mail: todorovt@vxcern.cern.ch 



TwoGAM implements Eq. ^ exactly for leptonic final states and for the direct, QPM, component 
of quark production, including all helicity states and the exact 2^4 kinematics. Single- and 
double-resolved components are added to the transverse-transverse scattering cross-section 
according to any parton distribution function selected from PDFLIB, and include a simple P^- 
suppression model. The resulting partonic states are hadronized by the Lund string model, 
as implemented in Jetset. The user can select from several different soft VDM scattering 
models. 



5.6 TWOGEN 



Author: A. Buijs, W.G.J. Langeveld, M.H. Lehto, D.J. Miller [|T|] 

E-mail: buijs@fys.ruu.nl 
Program size: 1800 lines 

Program location: http://www.fys.ruu.nl/ buijs/twogen/twogen.for 



TwoGEN samples the transverse-transverse luminosity function for real and virtual photons 
of the multiperipheral diagram. The events are then weighted with any user supplied cross- 
section (Jtt{i1 ^) and chosen using a simple "hit or miss" sampling. Thus, for example, 
the program can generate according to a chosen F2 or can produce resonances. Final state 
partons are fragmented using the Lund string model in Jetset. 
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5.7 Vermaseren 



Author: J. A.M. Vermaseren 

E-mail: t68@nikhef.nl 

The Vermaseren Monte Carlo has become the de facto standard calculation of e^e^ — >■ e"^e~ff 
via 77 collisions. It uses the subset of the exact 2 — >■ 4 matrix elements in which the electron 
and positron lines do not annihilate (i.e. formally they apply to the process e^/i~ — > e"'"/i~ff 
with equal electron and "muon" masses). 



6 Conclusions 



77 physics at LEP 2 is virtually impossible without event generators. Not only are they needed 
for modelling detector effects, as in almost all high energy physics analyses, but also because 
in 77 collisions the energy of the incoming photons is generally unknown. This means that we 
need to reconstruct the basic parameters of events solely from final-state properties. Except for 
the simplest final states consisting of only a few particles, detailed understanding is presently 
only possible through models implemented in event generators. 

During the course of this workshop, several general purpose QCD generators have been 
developed to better handle 77 and e7 interactions, enabling us to use experience gained from 
experiments with e"'"e~ , ep and pp collisions, where these programs have been used extensively. 
A lot more work needs to be done however. In contrast to the special purpose generators 
used so far in 77 physics, which in the case of high energy multi-particle production are less 
theoretically advanced, the general purpose ones have not yet been extensively confronted with 
available data. Such comparisons have already been started for 7p data from HERA with 
promising results, and it is important that this is also done with LEP 1 77 data in preparation 
for LEP 2. 
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